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Abstract 

We describe three different methods for generating quasi-exactly solvable potentials, 
for which a finite number of eigenstates are analytically known. The three methods are 
respectively based on (i) a polynomial ansatz for wave functions; (ii) point canonical trans- 
formations; (iii) super symmetric quantum mechanics. The methods are rather general and 
give considerably richer results than those available in the current literature. 
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In recent years, many authors have discussed examples of potentials for which a finite 
number of eigenstates are explicitly known, but not the whole spectrum. The first examples 
of such quasi-exactly solvable (QES) potentials were given by Razavyjjlj] and others[Q, ||. 
Subsequently, using a more general approach, many QES potentials were obtained by Tur- 
biner, Ushveridze and Shifman^, [5], [| and their connection with finite dimensional group 
representations was established [[?], |j. Also, Roy and VarshniQ gave an example of the 
generation of a new QES potential using supersymmetric quantum mechanics. 

In this paper, we develop and discuss three general methods for generating QES poten- 
tials. The first method is a systematic extension of the approach of Turbiner|5|, which has 
the distinctive feature of a polynomial ansatz of the same degree for all analytically known 
eigenstates, in addition to a conventional exponential factor for satisfying boundary con- 
ditions. Our extension provides a clearer motivation and consequently yields many three 
parameter QES potentials which are more general than those available in the literature. 
The method hinges on the choice of an initial variable z(x) and its use in making a specific 
ansatz for wave functions. 

The second method consists of making point canonical transformations (changes of 
dependent and independent variables) 0] on known QES potentials in order to obtain 
new ones fTlf . This approach is motivated from the observation that all known exactly 
solvable problems are inter-related via point canonical transformations. 

The third method discussed in this letter applies the techniques of supersymmetric quan- 
tum mechanics ]l2| to calculate the partner potentials of QES potentials with (n + 1) known 
eigenstates. As expected from the level degeneracy theorem for unbroken super symmetry, 
the supersymmetric partners are necessarily new quasi-exactly solvable potentials with n 
known eigenstates. 

Method 1: Polynomial Ansatz for Wave Functions: 

This method is a generalization of Turbiner's approach ||, [7| for generating QES poten- 
tials. We want to solve the Schrodinger equation (fi = 2m = 1) 

- tp"(x) + [V(x) - E] iJj(x) = , (1) 

along with the boundary condition that wave functions for bound states vanish as x ap- 
proaches the end points of the domain. The initial step in Method 1 for generating quasi- 
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exactly solvable potentials with (n + 1) analytically known eigenstates is to choose a non- 
singular monotonic variable z(x) and write the wave function ip in the form 

1>{x) = Pn (z, a ,ai,..., a n _x) e~ I y ^ dx , (2) 

where 

n 

Pn(z) = ^2cLjZ j (a n = 1) . (3) 
3=0 

is a polynomial of degree n in the variable z. [Note the special feature that a polynomial of 
the same degree n is present in all analytically known wave functions - however, different 
wave functions will correspond to different values of cij and consequently will have differing 
numbers of nodes]. The functional form of y(x) remains to be determined. Substituting for 
ip(x) into Eq. ([!]) yields 

V(x) -E = y 2 -y' + P " - 2yP ' n . (4) 

Pn 

In order that the potential have no poles in the physical domain, we require p n — 2yp' n to 
have p n as a factor. Also, in general, one expects that the energy will be a function of 
parameters in y(x), as well as the coefficients a,j in p n . However, we do not want these 
parameters to show up in the potential, as otherwise a dependence of the potential on 
energy would creep in. One way to prevent this from happening is to require that — ~ '' 



Pn 

be a polynomial of degree one. This condition restricts what the needed functional form for 
y[x) ought to be. More specifically, we choose y(x) by requiring that the quantity p' n — 2yp' n 
be a polynomial of degree (n + 1) in z. 

n+i 

p'n - typ'n = b i zJ ■ ( 5 ) 

This is equivalent to requiring 

n n n+1 

(z" - 2yz') ]T ajjzi- 1 + z' 2 £ a 3 j{j - l)z^ 2 = £ . (6) 

j=l 3=2 j=0 

For n = 1, the z' 2 term is absent, and it is sufficient to choose y(x) such that 

z" - 2yz' = q + qi z + q 2 z 2 . (7) 

Although the three constants gcb Qi, Q2 in y( x ) are arbitrary, their range of values is usually 
restricted by the requirement that the bound state wave function ip(x) vanish at the end 



points. For n > 1, the z' 2 term in Eq. (|6|) does contribute. Therefore, one needs the 
additional necessary constraint that the choice of initial variable z (x) be such that 

z' 2 = B + B x z + B 2 z 2 + B 3 z 3 . (8) 

This above constraint ensures that p'^ — 2yp' n is a polynomial of order of (n + 1). It also 
constrains the choices allowed for z(x), which now must have the form 

dz 



^—^^ = x (9) 

VB + B lZ + B 2 z 2 + B3Z 3 

The variables z = x ±2 , cos x, cosh x, e ±x , sech 2 x, ■ ■ •, etc. used in references ||, |||, [f| [Tj, are 
all special cases of the above Eq. @. Dividing Eq. (||) by Eq. (^) gives 

K - 2^. = + _ Qw _ i6w+i) + ^ n(z) (10) 

where the remainder polynomial is 

n-l 

^n( z ) = z3 fn{j) ; /n(i) = bj - aj_i6 n+ i - aj6 n + a,-a n _ifc n+ i . (11) 

j=0 

Clearly, the requirement of no poles in the potential V(x) implies 

/n(j)=0 , j = 0,l,...,n-l . (12) 

These n constraint equations can be solved to determine the constants ao, a%, . . . , a n —i. 
There are n solutions for any constant aj, since it satisfies an equation of degree n. The 
potential V{x) and energies E are given by 

V(x) = y 2 -y' + b n+1 z(x) ; E = -b n + a n _i6 n+1 . (13) 

The potential V(x) depends only on the parameters qo,Qi,Q2 appearing in y and b n+ \ = 
nq2 + n(n — 1)B 3 . Also, since b n = nq\ + (n— l)g2 a n-i + n(n — 1)B 2 + (n — l)(n— 2)i? 3 a n _i, 
a useful, alternative form for the energies is E = —nq\ — n(n — 1)B 2 + a n ^\[q 2 + 2(n — 1)^3]. 
Clearly, there are n energy eigenvalues corresponding to the n values of the constant a n _i. 

For concreteness, let us fully describe the important special case corresponding to the 
value n = 1. Here, we are constructing QES potentials with two analytically known eigen- 
states. The wave functions for these two states are both of the form given in Eqs. (|2|) and 

©■ 

ip(x) = Pl (z,a ) e-f yda , Pl (z) = a + z(x) . (14) 



The quantity y(x) is given by Eq. (0). For this case, the constraints given by Eqs. (12) 
read 

/i(0) = qo ~ a q 1 + a 2 q 2 = . (15) 

This quadratic equation for ao has two solutions. The QES potential is V{x) = y 2 — y' + 
q2z(x), and it has two analytically known eigenvalues 

E = a q 2 -qi = ~ ± -\J ' q 2 - 4:q q 2 , (16) 

with corresponding eigenfunctions given by Eq. (|l~4|). 

To illustrate the above-described polynomial ansatz method for generating QES poten- 
tials with (n+1) known eigenstates, we give three explicit examples with different choices of 
the initial variable z(x). These examples have been specifically chosen, since the potentials 
they generate have not been previously discussed in the literature. 
Example 1: z = x e ,n = 1 [2 known eigenstates]. 
We consider the polynomial p\ = z + ao- Eq. (0) then gives 



e 



1 1 / 

y = ~y — 2e \ q2X ^ ± + qiX + q ° x 

The QES potential with two known eigenstates is 

V 2 (x) = ^±^^-l[g 1 (e-2)-2g 2 ( e + l)^ + 2g (6-l)x- e ] 



n^ 1 _L „ „ _1_ r. ^-£+1 



1 

+ 4? 



2q mx x 2 -" + 2q iq2 x 2+e + (q 2 + 2q q 2 )x 2 + q 2 x 2 ~ 2e + q 2 x 2+2e . (17) 



The range of the potential is the half-line < x < oo [unless e = 1, in which case the range 
is the full line]. The known energy eigenvalues are given by Eq. (|l6|), and the eigenfunctions 
obtained using Eq. (113) are: 



tp = (x e + ao) exp 



1 — e t , 1 ( 12 _2+6 , gQ „2- e , 1 „ _2 



logx H H x z e + -gix' 



for e / ±2. (18) 



2 2eV2 + e 2-e 2 

The arbitrary constants are constrained by requiring that ifi satisfies boundary conditions. 
For instance, for any choice e > 2, one needs qo > 0, q 2 < 0. As an explicit example, 
consider e = 3, qo = l,q 2 = — l,q\ = 0. This gives ao = Tl ; eigenenergies i£ = ±1 and 
eigenfunctions 

(x 3 =F 1) r i/^ 5 -hi 

exp --(— + a; 1 . 

X DO 

Note that the polynomial x 3 + ao has no zero for the ground state, but one for the first 
excited state. 
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Example 2: z = coshx , [(n + 1) known eigenstates]. 

Since this choice of variable satisfies Eq. (Js|) , one can use it to generate QES potentials 
with an arbitrary number of known eigenstates. [The variable coshx also appears in Table 
4.1 of Ref. flj~2|| , and corresponds to the exactly solvable generalized Poschl- Teller potential 
V(x) = A 2 + (B 2 + A 2 + A)cosech 2 r - B(2A + 1) coth r cosech r with eigenenergies E n — 
A 2 - (A-n) 2 }. First we generate a quasi-exactly solvable potential with two eigenstates by 
taking the polynomial p±(x) = coshx + ao- This gives the following form for the function 
y(x): 

y(x) = a sinh x — /3cosechx — 7 coth x, (19) 
where a, /3, 7 are arbitrary constants. As before, ao obeys a quadratic equation 

2cm 2 , + (1 + 27) a - 2(a + 0) = 0. (20) 

Its two roots are given by 

-(1 + 2 7 ) ± V(l + 2 7 )^ + lto^+g) . 

a = • (21) 

These roots give rise to different number of nodes in the polynomial pi and hence in the 
wave functions ip: 

{ r\ 7-/3 / r \ 7+/3 
V'oc (coshx + a )e- QCOsh:E fcosh-J ( sinh 2j 

The requirement that i\) vanishes at x = 00 and x = implies that a > and 7 + (3 > 0. 
The potential is 

V2(x) = (j 2 — 7 + /3 2 ) coth 2 x — a (27 + 3) coshx + a 2 sinh 2 a? + ^(27 — 1) cothx cosech x. 

(22) 

It is independent of ao and singular at the origin, where it diverges as 1 ) _ xhe 

two lowest energy eigenvalues are given by 

E = ha/3 + (3 2 - 37 - 2aa - l) , (23) 



where, the solutions ao are given in Eq. (|21|). 

From Eq. (13), and using b n+ i = nq<2, we can get a potential with n + 1 states 



V n+ i = ^7 2 — 7 + j3 2 ^j coth 2 x— a (27 + 2n + 1) coshx+a 2 sinh 2 x+P(2j—l) cothx cosech x. 

(24) 



The energy eigenvalues are given by Eq. (|l^), where a n _i satisfies an algebraic equation of 
degree n + 1. In particular, for n = 2, the equation for a\ ia a cubic: 

-2a 2 a? - a(7 + 67)0? + 2(4a 2 + 4a/3 - 57 - 2 7 2 - 3)ai + 4(a + 2/3 + 2a 7 + 2/J7) = 0. 

Note that a special case of the QES potential of Eq. (|24"|) corresponding to /3 = 0, 7 = 0, 1 
was considered by Razavy|l|. The two choices of 7 give the even and odd eigenstates 
respectively. For these situations, the range of the potential is —00 < x < 00, since there 
is no singularity at x = 0. Also, in the limit a — * 0, this model goes to the generalized 
Poschl- Teller potential. 
Example 3: z = sech 2 x. 

This is another example where we can get a quasi-exactly solvable potential with more than 
two known states. The constants Bj in Eq. (||) are B$ = B\ = 0, Bi = 4, B^ = —4. The 
function y[x) is 

y = ^3 + cosech 2x + ( ^- — 1^ coth x + ^ cosh 2 a; coth x, 

and the potential is given by Eq. (^) with 6 n +i = ng2 — 4n(n — 1). A special case of 
this potential was discussed by Turbiner in Ref. ||. For n = 1, the two lowest energy 
eigenvalues are given by E = — q\ + q<i ao, and corresponding eigenfunctions are 

/ 2 \ \ ^ g2+6 ^ / \ -(12+11+10+ 2 ) 

ip = ( sech x + ao ) (cosh x) 4 (sinh x) 4 exp 



go cosh 2x 
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Method 2: Point Canonical Transformations: 



Point canonical transformations (PCT) are known [1C] to transform the Schrodinger 
equation for a potential V(x) into the Schrodinger equation for a new potential V(£). In 
this section, we employ PCT to generate new quasi-exactly solvable models starting from 
a known one. In a point canonical transformation, one changes from the independent 
variable x to the variable £, where x = /(£) and the wave function is transformed by: 
tp(ai]x) — > ^(ct-i'iO- The new potential and energy (denoted with tilde on the 

top) are then given by: 

1 / If" 2 f"'\ 

V-E = f 2 [V(a t ; /(f)) - E(a t )} + -I* - J —\ (25) 



The constant term on the right hand side represents the energy of the new potential. It 
should be emphasized that the above-described technique is quite general. We will illustrate 
it with a specific example. 
Example 4: 

As a starting point, we choose the new QES potential V^x) obtained in Eq. ( ]22| ) of Example 
2. The change of variables x = i£ converts hypergeometric functions into trigonometric 
functions. The new potential for < £ < tt is: 

V(£) = (7 2 - 7 + P 2 ) cot 2 £ + a (2 7 + 3) cos £ + a 2 sin 2 £ + 0(2j - 1) cot £cosec£. (26) 

The two known eigenvalues of this potential E are given by — E, with E given by Eq. (p3|). 
The corresponding wave functions are given by 

V> oc (cos£ + a )e- acos « 

The constraints 7 — /3>0,7 + /3>0 ensure that ip vanishes at the end points. 

Alternatively, if we make the change of variables x = coth _1 [isinh £], then the new 
potential is 

V = (3a-/3 + 2a7 + 2/3 7 )tanh£- (2a 2 -/3 2 + J B + 7-7 2 )tanh 2 £ 

- (3a + 207) tanh 3 £ + a 2 tanh 4 f . (27) 

This is a well defined QES potential over the range — cxd < £ < 00. The two lowest 
eigenvalues are E = —a 2 — E , where E is given by Eq. Q23| ) and the eigenfunctions are 

ifj cx (tanhz + a )e- atanh2 (tanhz) 7+/3 (sechz) 7 ~ /3 -^ . (28) 




Method 3: Supersymmetric Quantum Mechanics: 

If one has a potential V-{x) whose lowest eigenstate has energy E { ) = and wave 
function ipo(x), then its supersymmetric partner potential V+(x) is given by [^] 



w> = v -V>-*nKt) ■ (29> 
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The degeneracy theorem for unbroken supersymmetry states that V+ and V- have the same 
energy levels (except for the zero energy ground state) and simply related eigenfunctions: 

m = 4-» ; *t\ = 4- ri {i - fj) 4r ] ; a- = i. 2, . . . ( 3 o) 

Clearly, if one applies the above results to a QES potential with known eigenstates, 

one gets a new QES potential with n known eigenstates. For the simplest case of n = 1 and 
V-(x) = x 6 — 7x 2 + 2\/2, this procedure was discussed in Ref. ||. 

As a general illustration, we want to apply the supersymmetry approach to a QES 
potential obtained by Method 1 using a polynomial ansatz for the known wave functions. 
Since Eq. (||) gives 

Voocpoe"^ dx , 

the superpotential — V'o/Vo appearing in Eqs. (29) and ( |30| ) is just y(x) — p'o/po, where 
the parameters a,j appearing in the polynomial po correspond to the ground state. The 
supersymmetric partner potential is 

V + (x) = V_(x) + 2y'(x) - 2(P0P ° 2 " P ° 2) (31) 

Po 

and the lowest n eigenfunctions are 

oc ( & - &) p k e~ I , fc = 1, 2, . . . , n . (32) 

\Pk PoJ 

where pk is the polynomial p(z) with parameters aj corresponding to the /c th state of V-(x). 
Note that the same procedure of generating a supersymmetric partner potential can again 
be repeated by starting from the potential V+(x). Also, the standard methods of super- 
symmetric quantum mechanics ]l2| can be used to generate multi-parameter isospectral 
potential families of V-(x), which are of course new QES potentials. 
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